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Abstract 

The surface tension at the deconfinement transition of SU(3) is deter- 
mined with a parametrized version of the fixed point action of a renor- 
malization group transformation on lattices with temporal extent N T = 3 
and 4 and spatial extent N a /N T — 3 and 4. A considerable cut-off de- 
pendence can be seen in comparison with earlier determinations from tree 
level Symanzik and tadpole improved actions. 
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1 Introduction 



Lattice Monte Carlo simulations are at present the most effective tool to study 
SU(3) at finite temperature both in the gluon plasma phase, where perturbation 
theory is plagued by the well known infrared problem jjj], and at the critical 
temperature T c , i.e. in a highly non-perturbative regime. Lattice calculations 
of the bulk thermodynamic quantities in the deconfmcd phase of SU(3) with the 
Wilson action are affected, however, by a large cut-off dependence g], which is 
evident also in perturbation theory in the high temperature ideal gas limit Q . 
The continuum limit can be extrapolated with the Wilson action only from 
lattices with temporal extent N T not smaller than 6 j|] -0 Recent calculations j| 
U have revealed that such cut-off dependence can be drastically suppressed, from 
temperatures slightly above T c onwards, if tree level Symanzik Q or tadpole |j| 
improved actions are used in numerical simulations. In particular, with the 
tadpole improved action the free energy density shows no scaling violation from 
iV T = 3 to N T = 4 and is consistent with the continuum extrapolated from the 
Wilson action results. Moreover, the lattice determinations of the energy density 
minus three times the pressure from the tree level Symanzik improved action at 
N T = 4 are in agreement with the continuum. The latter result suggests that 
the cut-off effects beyond the tree level order are relatively unimportant in the 
high temperature phase of SU(3) jij. In other words, the tree level improvement 
alone can eliminate most of the cut-off dependence. 

Fixed point actions are lattice actions living in the space of couplings on 
the straight line which originates at the fixed point (FP) of a renormalization 
group (RG) transformation and leaves the critical surface along the direction of 
the coupling g 2 . They are classical perfect actions, i.e. their spectral properties 
are free of cut-off effects at the classical level. For this reason and in view 
of the previous considerations, FP actions are expected to yield a remarkable 
improvement in lattice studies of SU(3) thermodynamics. 

Given an arbitrary RG transformation, the FP action of any lattice config- 
uration is well defined and can be determined numerically by multigrid mini- 
mization 0, [|. For practical reasons, however, in Monte Carlo simulations only 
simple enough parametrizations of the FP action can be used. Being a potential 
source of cut-off effects, the parametrization of a FP action is a very delicate 
task. In SU(3) lattice gauge theory, the parametrizations proposed so far are of 
the form 

C i>l ^ 

2 In lattice simulations the temperature T and the volume V are determined by the lattice 
size N% X N T , N t < N a , through T = l/(N T a) and V = (A^a) 3 . 
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where the first sum is over all the closed paths C and Uc is the product of the 
link variables U^n) along the path C. The couplings Ci(C) are determined by 
a fit procedure on the numerically estimated "exact" values of the FP action of 
a representative set of Monte Carlo configurations |], [H]] . In Ref. jllj the arbi- 
trariness in the choice of the RG transformation has been exploited to define a 
very rotationally symmetric RG transformation, optimized in order that the re- 
lated FP action is very short-ranged and, therefore, easier to be parametrized. 
This RG transformation has been called for historical reasons "type III" RG 
transformation. The related FP action has been parametrized with the pla- 
quette and the twisted perimeter-6 loop (see Fig. [I]), with the couplings Ci(C), 
i = 1, .... 4, given in Table |. 




plaquette twisted perimeter-6 

Figure 1: Loops considered in the type III parametrized FP action. 

In Rcf. |nj the type III parametrized FP action has been used to determine 
from the lattice the free energy density of SU(3) at T/T c = 4/3, 3/2, 2. Sim- 
ulations were performed on lattices as small as 8 3 x 2 and 12 3 x 3. Results at 
N T = 3 are in agreement with the continuum as extrapolated with the Wilson 
action from lattices with N T — 6 and 8 § and with the results from the tadpole 
improved action at N T = 3 and 4 (for a complete review, see also Q ). 



Table 1 : Couplings of the type III parametrized FP action [[To| . 





Cl 


C2 


C3 


c 4 


plaquette 


0.4822 


0.2288 


-0.1248 


0.0228 


twisted perimeter-6 


0.0647 


-0.0224 


0.0030 


0.0035 



In this paper, the same parametrized FP action has been used to determine 
the surface tension of SU(3) on lattices with temporal extent N T — 3 and 4 
and spatial extent N a /N T = 3 and 4. Like the latent heat, the surface tension 
is a physical quantity characteristic of the first order deconfinement transition 
of SU(3)J^] As a consequence, it is affected by the properties of both the low 
and the high temperature phase. Since in the latter phase the use of the type 

3 If the high temperature phase transition of the full QCD were first order, the surface 
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Ill parametrized FP action allows an impressive reduction of the cut-off depen- 
dence of the thermodynamic quantities, it is interesting to check if the same 
occurs for a quantity typical of the critical region. 

The results for the surface tension obtained in this work have been then com- 
pared with earlier determinations from the Wilson action jl3| , tree level Syman- 
zik and tadpole improved actions [fuj|. 



2 Monte Carlo results for the surface tension 



At the critical temperature of the first order deconfinement transition of SU(3), 
there can be mixed states where the confined and the deconfincd phases coexist, 
separated by an interface. These mixed states have an additional free energy 
F = a A (a is the surface tension, A the area of the interface), which indicates 
that they are less probable than pure states of one or the other phase. The 
frequency distribution of any order parameter £1 at the transition has a typi- 
cal double-peak structure, where the two peaks correspond to the pure phase 
configurations, while the region in-between corresponds to configurations con- 
taining an interface. The peaks become more pronounced when the volume is 
increased. The probability distribution of fl is given by 

P(Q) = Cl e~^ T e~^' n ^/ d " + c^-hvlT^-iSi^fldl 

+ c m c~ {hVl+hV ' 2+ ' 2 ' TA)/T , (2) 
Cl ocy 1/2 , d t ctV- 1/2 , i = 1,2 , 

where /i, /2 are the free energy densities of the two pure phases and Vi, V% 
are their volumes. The factor 2 in the free energy of the interface appears since 
at finite volume with periodic boundary conditions two interfaces are needed 
to separate two volumes. In the previous formula, the dependence on Q enters 
through the relation £1 = (Viili + Vi^l^jV ■ At T c and at infinite volume, the 
free energy densities in both phases are identical, /i = f%. The leading volume 
dependence of the surface tension is determined by 



a \ 1 /A r \ 2 ] / P x 



T-.-i ) •) \ v I m p7l p72 I J ( 3 ) 

± c J V V v <r / y r max,l r max,2 / 

where P mm is the minimum of the distribution P{£1), -P m ax,i and -P max .2 are 
the two maxima corresponding to the values f2i and f2 2 of £1 in the pure states 
of the two phases at infinite volume. The weights 71 and 72 are determined 
imposing (ft) = jiQi + 72^2, with 71 + 72 = 1- 

tension would play an important role in the quark-gluon plasma formation in heavy-ion ex- 
periments and in the nucleation of hadronic matter in the early universe |l3[ . 



4 



Table 2: Parameters of the runs and results for j3 c and xl/N%. Errors have 
been estimated by the jackknife method. 



lattice 


#/? 


# iterations 


P reweighting 


I3 C 




12 3 x 3 


3 


645099 


3.58995 


3.58982(9) 


2.958(20) x IQ-' A 


12 3 x 4 


1 


171209 


3.69915 


3.7007(3) 


1.013(16) x 10~ 2 


16 3 x 4 


3 


135662 


3.70025 


3.7009(5) 


8.87(20) x 10~ 3 



A convenient choice for the order parameter in SU(3) is the absolute value 
of the Polyakov loop L = 1/./V 3 J^n TrJln^i U^(ri, n 4 ). Numerical simulations 
have been performed on three lattices at values of f3 close to criticality. A 
summary of the parameters of the runs is given in Table |^. Here one iteration 
means the combination of a 20-hit Metropolis and 4 over- relaxation updatings; 
the runs have been done on DEC-alpha machines. In all the numerical simu- 
lations frequent flips of the order parameter \L\ could be observed between the 
two phases, thus indicating a proper sampling of the thermal ensemble at criti- 
cality. The occurrence of frequent jumps was expected in view of the relatively 
low N a /N T ratio. The critical couplings have been determined through the lo- 
cation of the peak in the Polyakov loop susceptibility xl — A 3 ((|iy| 2 ) — (|L|) 2 ). 



They are in agreement within errors with the results of Ref. 1 10 where the so 
called Columbia definition was adopted^] The Polyakov loop distribution on the 
12 3 x 3 lattice and on the two lattices with N T — 4 is presented in Fig. ||. Data at 
different ft on the same lattice have been interpolated by Ferrenberg-Swendsen 
reweighting in order to make the peaks in the Polyakov loop distribution have 
the same height. In order to determine P m i n , -P m ax,i and P max ,2, the statistical 
fluctuations in the vicinity of the extrema of P(\L\) have been smoothed out by 
a polynomial fit. 

In Table |3| the results of this work for ct/T 3 are compared with other deter- 
minations from the Wilson action |i~3] ] , the tree level Symanzik improved action 
(plaquette +(1x2) loop) and the tadpole improved action (with the same 
loops) p3 . The infinite volume extrapolations have been done according to the 
ansatz [|13l 



T?) v \Tr ) \N n 



(4) 



The result on the 12 3 x 3 lattice and the infinite volume extrapolation at N T = 4 
for the tadpole improved action indicate no cut-off dependence and a contin- 
uum value of er/T 3 equal to 0.0155(16) (corresponding to a ~ 7 MeV/fm 2 ).^ 

4 This comparison is important since in Rcf. the physical scale in the determination of 
the free energy density was set using the critical couplings as obtained in Ref. [juj. 

5 Here and in the following it is assumed that the infinite volume extrapolation at N T = 3 
does not differ too much from the result on the 12 3 X 3 lattice. 
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Table 3: cr/T^ for various lattice actions on several lattices. Column 7 gives 
the results for 71 from the type III parametrized FP action. Errors have been 
estimated by the jackknife method. 







Wilson |3| 


tree level JLJ 


tadpole JTj 


type III FP 




volume 








7i 


3 


12 3 




0.0234(24) 


0.0158(11) 


0.0307(8) 


0.420(10) 


4 


12 2 x 24 


0.0241(27) 










4 
4 
4 
4 
4 


24 2 x 36 
12 3 
16 3 
24 3 
32 3 


0.0300(16) 


0.0148(16) 
0.0136(25) 
0.0116(23) 


0.0147(14) 
0.0119(21) 
0.0125(17) 


0.0196(11) 
0.0180(21) 


0.409(21) 
0.441(23) 


4 


00 


0.0295(21) 


0.0152(26) 


0.0152(20) 


0.026(5) 




6 
G 
6 


20 3 
24 3 
36 2 x 48 


0.0123(28) 
0.0143(22) 
0.0164(26) 
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00 


0.0218(33) 











Moreover, this continuum value is in agreement with the infinite volume ex- 
trapolation at iV r = 4 from the tree level Symanzik improved action and with 
the extrapolation in l/N^ to N T = 00 of the results with the Wilson action at 
N T = 4 and 6. 

In the case of the type III parametrized FP action, the result on the 12 3 x 3 
lattice and the infinite volume extrapolation at N T — 4 are consistent with scal- 
ing, but at a value considerably larger than the continuum determined from 
the tadpole improved action. If one observes, however, that the infinite volume 
extrapolation at N T = 4 is probably not reliable for the smallness of the lattices 
involved and assumes that the true extrapolated value cannot be too different 
from the result on the 16 3 x 4 lattice, then it must be concluded that there is a 
large scaling violation. Although the result on the 16 3 x 4 lattice is not incon- 
sistent with the corresponding determinations from the tree level Symanzik and 
the tadpole actions, the result on the 12 3 x 3 lattice differs from the continuum 
more than in the case of the tree level Symanzik action. 

This scenario would suggest that the type III parametrized FP action does 
not improve at T = T c , at least as far as the surface tension is concerned, 
although it clearly does, from temperatures slightly above T c onwards, in the 
case of the free energy density [|ll] . 

Before accepting this conclusion, possible sources of error have been investi- 
gated. The procedure followed to find the results of this work has been verified 
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by re-calculating some old results for the surface tension at iV T = 2 obtained 
with the Wilson action Jl5[ . 

The use of the naive Polyakov loop as order parameter, instead of the cor- 
responding classically perfect operator, could represent in principle a source of 
lattice artifacts. This possibility has to be excluded, however, for several rea- 
sons: first of all, the type III parametrized FP action was optimized in Ref. jlO| 
in order to dump the violation of rotational symmetry in the correlators of 
naive Polyakov loops and, consequently, in the static qq potential. The tests 
performed in Ref. provided impressively good results, in comparison with 
previous parametrizations of the FP action. Moreover, the surface tension de- 
pends only on the height of the minimum and of the maxima of the distribution 
P(\L\) 1 while possible lattice artifacts in the corresponding values of \L\ should 
play a minor role. This statement is supported by the fact that the nice re- 
sults for the surface tension from the tadpole improved action of Ref. jl4| were 
obtained with the naive Polyakov loop. Finally, it has been explicitly checked 
in the present work that using the FP action density itself as order parameter 
instead of the absolute value of the Polyakov loop does not change the scenario, 
although the statistical fluctuations are larger. 

The most natural explanation of the large observed cut-off dependence is 
that it is induced by the parametrization of the FP action, which perhaps fails 
to represent the exact FP action on configurations typical of the critical region. 
By comparing the results from the type III parametrized FP action with those 
of Ref. [Q, one could also argue that planar (although less local) loops are 
preferable in lattice actions for simulations in presence of surfaces. The fact that 
the type III parametrized FP action does not satisfy the Symanzik condition at 
the tree level |l(| could play a role in this respect. 

It must be stressed, however, that the surface tension is a difficult quantity 
to determine for the high statistics required and for the complicated volume 
dependence. Further checks are necessary to confirm the previous conclusion. 
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Figure 2: Polyakov loop distribution on the 12 3 x 3 lattice and on the 12 3 x 4 
and 16 3 x 4 lattices. Data on the different lattices have been normalized to the 
corresponding total statistics. 
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